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A  new  dimensionless  relationship  for  analysis  of  ballistic  penetration  data  is  derived  and  applied  to 
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per  unit  mass.  Simultaneous  consideration  of  results  of  the  dimensional  analysis  and  material  properties 
suggest  that  the  length  scale,  which  is  related  to  decreasing  penetration  depth  with  increasing  tile 
thickness,  correlates  with  the  ratio  of  surface  energy  to  elastic  modulus.  The  energy  per  unit  mass,  which 
is  linked  to  the  relationship  between  penetration  depth  and  penetrator  velocity,  correlates  with  dynamic 
shear  strength  of  failed  ceramic  reported  from  plate  impact  experiments,  divided  by  mass  density.  The 
dimensional  analysis  provides  a  structured  framework  under  which  future  multiscale  simulations  and 
validation  experiments  can  be  organized  and  compared. 
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1.  Introduction 

Ceramic  materials  comprise  important  components  of  many 
modern  protection  systems  for  vehicles  and  personnel,  i.e.,  vehic¬ 
ular  armor  and  body  armor.  Experiments  and  numerical  simula¬ 
tions  of  ballistic  performance  of  armor  ceramics  have  considered 
various  projectile  and  target  geometries,  materials,  and  velocity 
regimes,  as  summarized  in  Ref.  [1,2],  for  example.  Popular  experi¬ 
ments  include  those  that  measure  limit  velocities  for  complete 
perforation  of  monolithic  or  layered  structures  (i.e.,  V50),  those 
intended  to  monitor  dwell  in  confined  cylindrical  targets,  and  those 
that  measure  penetration  depth  into  a  backing  material,  typically 
metallic,  placed  behind  ceramic  tile(s).  Properties  often  deemed 
favorable  include  high  hardness,  high  elastic  stiffness,  high 
strengths  (static/dynamic  compressive,  shear,  and  bending),  and 
low  density  relative  to  armor  steels;  drawbacks  include  brittleness 
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and  high  cost.  However,  consensus  regarding  relative  importance  of 
ceramic  material  properties  affecting  results  of  experiments 
designed  to  inform  protection  applications  remains  elusive.  Re¬ 
lationships  among  macroscopic  properties  and  microstructure 
features  are  difficult  to  isolate  and  hence  often  unclear. 

Efforts  towards  relating  measurable  properties  to  ballistic  per¬ 
formance  of  brittle  materials  (e.g.,  ceramics  and  glass)  have  been 
undertaken  since  at  least  the  1960s.  Comprehensive  review  articles 
on  the  subject  in  the  context  of  armor  ceramics  include  [1,3,4];  a 
few  other  notable  works  are  mentioned  here.  In  Ref.  [5],  correlation 
of  performance  with  an  average  of  the  static  and  dynamic 
compressive  strengths  of  the  ceramic,  normalized  by  mass  density, 
was  reported.  In  Ref.  [6],  it  was  suggested  that  ceramic  hardness, 
fracture  toughness  (or  surface  energy),  and  penetration  resistance 
may  be  related.  In  contrast,  recent  experiments  [7]  on  various  sil¬ 
icon  carbides  found  depth  of  penetration  to  depend  on  hardness 
but  not  toughness.  In  Ref.  [8],  dynamic  shear  strength  was  inferred 
as  an  important  property  in  ballistic  penetration  experiments; 
subsequent  computational  results  [9]  and  data  analysis  [10]  have 
further  demonstrated  important  effects  of  friction  and  shear 
strength  of  the  fractured  or  comminuted  ceramic  material  when 
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impacted  in  the  hypervelocity  regime.  In  Ref.  [11],  flaw  size,  elas¬ 
ticity,  plasticity,  and  fracture  properties  enter  an  analytical  model 
for  dwell  based  on  a  dimensionless  ductility  parameter  [12]. 
Importance  of  both  hardness  and  ductility  (including  inelastic 
deformation  mediated  by  dislocations  and/or  micro-cracks)  on 
transition  velocities  for  dwell  in  confined  thick  ceramic  targets  was 
reported  in  Ref.  [13]. 

Principles  of  dimensional  analysis  were  formally  established 
upon  proof  of  Buckingham's  pi  theorem  of  the  early  20th  century 
[14].  Dimensional  analysis  enables  a  reduction  in  the  number  of 
independent  variables  that  must  be  considered  in  a  functional 
relationship  among  physical  quantities  observed  in  experiments. 
For  example,  if  fundamental  units  of  mass,  length,  and  time  are  all 
involved  independently,  a  reduction  in  the  number  of  independent 
variables  by  three  is  possible.  Detailed  discussion  and  example 
applications  in  solid  and  fluid  mechanics  can  be  found  in  Refs. 
[15-17],  with  the  latter  [17]  formalizing  regimes  in  which  self¬ 
similarity  may  describe  a  particular  problem.  Dimensional  anal¬ 
ysis  may  enable  relating  data  collected  from  laboratory  experi¬ 
ments  conducted  at  one  length/time  scale  to  data  collected  from 
observations  of  physical  phenomena  at  scales  exceeding  pragmatic 
reach  of  experimental  facilities.  It  is  noted,  however,  that  the  choice 
of  dimensionless  parameters  is  not  unique,  and  physical  insight  is 
often  obtained  only  from  particular  choices  guided  by  logic  and 
some  additional  knowledge,  beyond  the  numerical  data,  of  the 
problem  under  consideration. 

The  only  previously  published  attempt  incorporating  formal 
dimensional  analysis  to  relate  ceramic  material  properties  to  bal¬ 
listic  performance  in  depth  of  penetration  experiments  known  to 
the  present  author  is  reported  in  Ref.  [18],  wherein  no  strong  cor¬ 
relations  were  discovered.  However,  this  prior  effort  considered 
only  a  small  scope  of  data  available  in  years  up  to  1992  and  focused 
primarily  on  static  (but  not  dynamic)  material  properties.  Perhaps 
further  complicating  (if  not  confusing)  the  analysis  are  simulta¬ 
neous  consideration  of  multiple  performance  metrics  (e.g.,  pene¬ 
tration  depth  normalized  in  various  ways)  that  lead  to  different 
rank  orderings  of  material  superiority  for  the  same  set  of  raw  data. 
Complicating  the  general  problem  of  analysis  of  ballistic  penetra¬ 
tion  data  on  armor  ceramics  is  the  large  number  of  factors 
[2,19-23]  that  may  affect  performance  in  addition  to  impact  ve¬ 
locity,  ceramic  material,  and  total  ceramic  target  thickness, 
including  but  not  necessarily  limited  to  the  following:  penetrator 
geometry,  penetrator  material,  backing  material,  possible  lateral 
conflnement,  possible  face-plate  or  buffer,  layers  of  thin  versus 
thick  tiles  with  possible  binder  material,  lateral  dimensions  of  tiles, 
and  variability  in  ceramic  material  purity  and  initial  defect  content. 
Dimensional  analysis  has  been  applied  towards  penetration  of 
mononlithic  metallic  targets  by  long  rods  with  mixed  success 
[2,24-26].  Recently,  similitude  analysis  has  been  used  to  describe 
penetration  of  steel  targets  by  tungsten  alloy  and  hardened  steel 
projectiles  [27],  including  development  of  equations  for  predicting 
ballistic  limit  thickness  and  limit  velocity,  as  well  as  residual  pro¬ 
jectile  velocity  and  residual  projectile  length  for  overmatched 
targets. 

The  present  paper  offers  the  following  new  technical  contribu¬ 
tions.  Recognizing  that  penetration  depth  versus  ceramic  target 
thickness  data  tend  to  fall  on  lines  of  constant  slope  regardless  of 
impact  velocity,  an  additively  separable  relationship  for  normalized 
depth  of  penetration  is  proposed.  This  relationship  contains  two 
material  parameters:  a  length  scale  ^  and  an  energy  per  unit  mass 
S.  The  relationship  is  fit  to  penetration  data  from  six  experimental 
studies  reported  in  the  open  literature  [5,23,28-31]  on  six  poly¬ 
crystalline  armor  ceramics:  low-  and  high-purity  alumina 
(respectively  labeled  AI2OT  and  A^O^),  aluminum  nitride  (AIN), 
boron  carbide  (B4C),  silicon  carbide  (SiC),  and  titanium  diboride 


(TiB2).  In  other  words,  parameters  and  S  are  determined  for  each 
set  of  experimental  data;  different  values  are  obtained  among  data 
sets  because  of  different  test  geometries,  backing  materials,  etc. 
These  two  parameters  are  correlated  with  ceramic  material  prop¬ 
erties  upon  consideration  of  physics  associated  with  each  param¬ 
eter  via  the  way  it  enters  the  dimensionless  equation.  Results  offer 
new  insight  into  performance  evaluation  of  armor  ceramics  in 
ballistic  impact. 

The  main  content  of  this  paper  is  organized  as  follows.  General 
concepts  from  dimensional  analysis  applicable  to  problems  in 
penetration  mechanics  are  considered  in  §2.  Particular  applications 
of  these  concepts  towards  ballistic  penetration  data  on  ceramics  are 
addressed  in  §3.  Relationships  among  parameters  entering  the 
dimensionless  equation(s)  and  properties  of  ceramic  materials  are 
postulated  in  §4,  with  corresponding  discussion.  Conclusions  follow 
in  §5.  Regarding  notation,  mathematical  variables  are  written  in 
italic  or  Greek  font,  with  a  definition  of  each  variable  given  soon 
after  it  first  appears  in  the  text. 

2.  Dimensional  analysis:  general 

In  this  paper,  Buckingham's  pi  theorem  is  applied  towards 
dimensional  analysis  of  physical  systems  involving  ballistic  pene¬ 
tration.  Firstly,  Buckingham's  pi  theorem  applied  to  any  system  can 
be  stated  as  follows  [14-16].  If  an  equation  involving  n  variables  is 
dimensionally  homogeneous,  it  can  be  reduced  to  an  equation 
among  n-k  independent  dimensionless  products,  with  k  being  the 
number  of  independent  reference  dimensions  needed  to  charac¬ 
terize  all  of  the  variables.  In  dimensional  (rather  than  dimension¬ 
less)  form,  let 

y  =/(X2.X3,....Xn),  (2.1) 

where  y  is  the  dependent  variable  and  X2,...,  Xn  are  independent 
quantities,  some  of  which  may  be  independent  variables  and  others 
constants  in  a  given  problem.  Equation  (2.1)  can  be  transformed 
into  dimensionless  form  as 

Til  =(/)(il2,7l3,...,iI„_/,),  (2.2) 

where  J7i  is  the  dimensionless  version  of  y  and  function  0  depends 
on  n  -  /<  -  1  other  dimensionless  products  (i.e.,  pi  terms)  con¬ 
structed  from  the  original  set  {y,  X2,  X3,...,  Xn}.  Pragmatically, 
normalization  of  all  variables  entering  tp  should  be  performed  using 
only  the  set  {X2,  X3, . . . ,  Xn}  so  that  the  independent  variable  y  appears 
only  once,  on  the  left  side  of  the  equation. 

The  ballistic  penetration  problem  to  which  Buckingham's  the¬ 
orem  is  applied  herein  is  shown  in  Fig.  1(a).  The  experimental  set¬ 
up,  as  introduced  for  example  in  Ref.  [5],  consists  of  a  projectile  of 
initial  length  Lq  impacting  a  ceramic  target  of  total  thickness  h  at 
velocity  Vq.  with  the  ceramic  target  fronting  a  thick  (effectively 
semi-inflnite)  metallic  backing.  The  ceramic  target  may  consist  of 
one  or  more  layered  tiles.  The  present  analysis  is  restricted  to 
normal  impact,  i.e.,  effectively  null  obliquity.  Performance  of  the 
ceramic  is  measured  by  residual  penetration  depth  P,  with  pene¬ 
tration  resistance  decreasing  with  increasing  P. 

In  dimensional  form,  penetration  depth  can  be  written  as  the 
following  function  of  impact  velocity,  geometric  variables  {g},  and 
material  property  variables  {m},  letting  subscripts  0,  T,  B  denote 
penetrator,  ceramic  tile,  and  backing  metal: 

P  =  P{Vo,  {S}o,  {m}o.  telr.  {mil.  {sIb-  {^Ib)-  (23) 

Possible  lateral  conflnement  is  included  in  the  sets  of  geometric 
variables  for  the  tile  and  backing.  In  subsequent  analysis,  applica- 
tion/fltting  of  (2.3)  is  restricted  to  experimental  data  sets  for  which 
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ceramic  tile(s) 


Fig.  1.  Ballistic  experiment  (a)  and  typical  fits  to  residual  penetration  data  (b). 


the  backing  and  tile  geometry  remain  fixed  and  the  penetrator 
material  remains  fixed.  Further  letting  {g}o  ^  Lq,  {g}r  ^  h,  (2.3) 
reduces  under  these  assumptions  to 

P  =  P(Vo,Lo,h,{m}r,mB)-  (2.4) 

A  reduced  form  of  (2.4)  is  needed  in  the  context  of  dimensional 
analysis.  This  reduced  equation  should  satisfy  the  following  re¬ 
quirements:  (i)  it  should  be  a  dimensionless  function  of  dimen¬ 
sionless  quantities,  (ii)  it  should  yield  the  penetration  depth  into 
the  bare  backing  material  as  h  ^  0,  and  (iii)  it  should  satisfy 
observed  physics  of  the  problem,  notably  a  linearly  decreasing 
penetration  depth  with  increasing  tile  thickness,  with  the  slope  of  P 
versus  h  constant  among  test  data  for  different  impact  velocities 
with  all  other  quantities  (i.e.,  material  properties  and  geometry) 
held  fixed.  The  latter  point  is  demonstrated  visually  in  Fig.  1(b),  and 
is  characteristic  of  ceramic  penetration  data  as  noted  in  Ref.  [2]  and 
demonstrated  later  in  results  shown  in  §3.  In  order  to  satisfy  these 
requirements  with  minimal  complexity,  three  parameters  are 
introduced:  a  reference  velocity  Vr  that  is  used  to  normalize  the 
impact  velocity  in  the  penetration  depth  relation  for  the  bare 
backing  material,  a  material  length  scale  ^  that  depends  on  the 
ceramic  tile  material,  and  an  energy  per  unit  mass  ^  that  also  de¬ 
pends  on  the  ceramic  tile  material.  In  other  words,  {m}T 
and  {mjB  Vr.  In  dimensional  form. 


P/Lo  =  4>  =  vUi/se)  +  (2.8) 

where  cp  is  linear  in  h,  (/)(0)  =  0,  and  \piVolVR,0)  =  Application  of 
(2.8)  to  data  in  §3  demonstrate  that  a  linear  penetration  depth 
versus  impact  velocity  relation  is  sufficient  to  describe  the  bare 
backing  metal  over  velocity  regimes  of  interest,  as  will  be  shown 
explicitly  later.  The  final  form  of  (2.6)  used  herein  thus  becomes 


P_ 

4 


h  Vo 


(2.9) 


Flere,  Cq  is  a  dimensionless  constant  entering  the  aforementioned 
linear  fit  to  penetration  data  for  the  bare  backing.  Penetration 
depth  decreases  linearly  with  increasing  tile  thickness  and  in¬ 
creases  with  increasing  Increases  in  ceramic  property  S  linearly 
reduce  the  velocity  dependent  component  of  residual  penetration. 
Consider  (2.8)  and  (2.9)  applied  to  a  data  set  for  a  fixed  tile  material, 
in  which  only  Vo  and  h  vary  among  data  points.  The  penetration 
depth  equation  is  an  example  of  a  separable  function  of  rank  p  =  2 
of  the  general  form  [32] 


mVo)  =  J2vM4'i{Vo), 


1=1 


(2.10) 


P  =  P{Vo.Lo,VR.h,^,^).  (2.5) 

Each  application  of  (2.5)  is  further  restricted  to  data  sets  for 
which  the  penetrator  geometry  is  fixed,  such  that  Lq  =  constant  is 
used  only  to  normalize  P  as  is  conventional  in  analysis  of  ballistic 
data  [2,24,25],  leading  to  the  dimensionless  form  [requirement  (i)] 

P/Lo  =  4>(Vo/Vr,  h/se,  W I Vl) .  (2.6) 

In  the  context  of  Buckingham's  theorem,  two  independent  di¬ 
mensions  (length  and  velocity)  enter  the  variables  considered  in 
(2.5),  so  k  =  2.  Since  Lq  is  fixed  and  has  units  redundant  with  h,  the 
total  number  of  pi  terms  is  n  =  6,  and  the  number  of  dependent 
dimensionless  terms  is  thus  n  -  /<  -  1  =  3,  which  is  fully  consistent 
with  the  number  of  dependent  variables  on  the  right  side  of 
dimensionless  Equation  (2.6).  Requirement  (ii)  above  implies 


(P/Lo)b  =  HVo/Vr,  0, 0)  =  0b(Vo/V^),  (2.7) 

where  05  is  the  dimensionless  penetration  depth  into  the  pure 
backing  metal.  Requirement  (iii)  suggests  a  separable  equation  of 
the  following  form: 


where  in  particular  here 


(Pl=(p  =  -h/^,  (P2  =  ^]  V^l=l,  ^2=^ 


Vr 


=  C0  +  T7^  1-777  • 


Vr 


(2.11) 


The  following  properties  then  apply  for  partial  derivatives  of  0: 


90  dcpi  dcp  90  90J  90 

dh  dh  *  dh^  9Vo  dVg 


(2.12) 


For  the  linear  functions  (p{h)  and  0(Vo)  in  (2.11),  Equation  (2.12) 
can  be  solved  for  dimensional  ceramic  material  parameters  ^  and 


Equations  in  (2.13)  can  be  derived  directly  from  (2.9)  for  the 
present  case  wherein  P  =  P(h,Vo)  with  other  factors  held  constant. 
Flowever,  (2.10)-(2.12)  imply  how  the  approach  could  be  extended 
to  other  more  general  kinds  of  separable  functions.  Note  also  that 
0  ^  in  (2.12)  unless  0  =  0. 
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3.  Dimensional  analysis:  application  to  test  data 

For  the  backing  metal  alone  (i.e.,  no  tile),  (2.9)  reduces  to 

P/Lo  =  Co  +  Vo/l/^.  (3.1) 

A  linear  regression  with  fit  quality  denoted  conventionally  by 
is  applied  to  determine  intercept  Cq  and  reference  velocity  (inverse 
slope)  Vr  for  each  data  set  commensurate  with  experiments  listed 
in  Table  1.  Resulting  values  are  shown  in  Table  3.  Sufficiency  of  the 
linear  fits  is  demonstrated,  for  example,  by  representative  results 
shown  in  Fig.  2.  The  quality  of  fits  for  which  only  two  experimental 
data  points  are  reported  [28]  is  indeterminate. 

Experimental  configurations  for  which  the  final  dimensionless 
form  of  residual  penetration  depth  equation  in  (2.9)  is  applied  are 
listed  in  Table  1.  These  consider  various  projectile  types  of  lengths 
ranging  from  12.7  to  72.5  mm  and  impact  velocities  ranging  from 
0.6  to  3.0  km/s.  Backing  metals  consist  of  monolithic  thick  plates  of 
various  steels  and  aluminum  alloys.  Ceramic  tiles  include  six 
different  effectively  isotropic  polycrystalline  materials:  low  purity 
alumina  (e.g.,  AD-85),  high  purity  alumina  (e.g.,  AD-995), 
aluminum  nitride,  boron  carbide,  silicon  carbide,  and  titanium 
diboride.  Total  tile  thickness  ranges  from  1.3  to  100  mm.  In  some 
cases  [28,30],  experiments  considered  stacks  of  multiple  thinner 
tiles  of  the  same  total  thickness  h  as  one  or  fewer  thicker  tiles;  in 
such  cases,  the  number  of  tiles  did  not  appear  to  affect  the  ballistic 
penetration  depth,  in  contrast  to  results  reported  in  one  different 
study  [22]. 

Physical  properties  of  interest  are  listed  in  Table  2  with  sup¬ 
plementary  references.  Initial  mass  density  p,  elastic  (Young's) 
modulus  £,  Poisson's  ratio  v,  fracture  toughness  Kc,  compressive 
strength  jq  bending  strength  ap,  and  Vicker's  hardness  Hy  are  ex¬ 
amples  of  static  material  properties.  The  Flugoniot  Elastic  Limit 
(FIEL)  an  and  dynamic  shear  strength  r  of  the  shocked  ceramic,  the 
latter  defined  in  Ref.  [10]  from  the  intersection  of  the  elastic  line 
with  the  failed  strength  curve  of  the  shocked  material,  are  exam¬ 
ples  of  dynamic  material  properties.  Lattice  parameters  are  a  and  c. 
Although  values  of  properties  shown  in  Table  2  are  deemed  most 
physically  representative  of  the  ceramics  tested  experimentally, 
complete  property  sets  require  consideration  of  characterization 
data  from  multiple  sources,  and  variations  among  properties  for  a 
given  material  are  common  as  evidenced  by  ranges  of  static  prop¬ 
erty  values  reported  in  Ref.  [18],  for  example.  All  of  these  ceramics 
undergo  fracture  (transgranular  and/or  intergranular)  when  sub¬ 
jected  to  dynamic/shock  loading  of  sufficient  magnitude.  Other 
known  inelastic  deformation  mechanisms  are  listed  for  each  ma¬ 
terial  in  the  second  column  from  the  right.  In  particular,  slip  and 
twinning  refer  to  dislocation  glide  and  mechanical/deformation 
twinning  [33],  respectively.  The  phase  change  in  AIN  is  a  solid- 
-solid  transformation  from  hexagonal  to  cubic  crystal  structure 
that  occurs  at  pressures  around  20  GPa  (more  specifically,  wurtzite 
to  rocksalt  [34]).  Amorphization  in  B4C  is  a  stress-induced  change 
from  trigonal  (i.e.,  rhombohedral)  crystal  structure  to  a  non¬ 
crystalline  solid  phase  (i.e.,  glass)  thought  to  occur  at  ballistic 
impact  stresses  in  excess  of  23  GPa  [35].  While  fracture  and  other 


inelastic  deformation  mechanisms  in  Table  2  presumably  affect 
dynamic  properties  (e.g.,  r),  reliable  quantitative  relationships  be¬ 
tween  thresholds  for  microscopic  mechanism  activation  and  values 
of  macroscopic  properties  remain  elusive  for  most  polycrystalline 
armor  ceramics. 

Equation  (2.9)  is  applied  separately  to  penetration  data  for  each 
ceramic  tile  material  considered  in  each  experimental  data  set 
listed  in  Table  1.  Specifically,  material  parameters  ^  and  S  are  fit  for 
each  such  application  as  follows.  A  plot  of  P  versus  h  provides  an 
estimate  of  0P/0h,  which  is  found  to  be  fairly  constant  for  each 
impact  velocity  Vo-  A  plot  of  P  versus  Vo  at  fixed  values  of  h  gives  an 
overall  estimate  of  0P/0VO,  which  is  found  to  be  fairly  constant  and 
not  vary  too  strongly  with  thickness  h.  Then,  Equations  in  (2.13)  are 
used  to  determine  an  initial  guess  of  each  of  3  and  S.  The  two 
parameters  c^and  s  are  then  adjusted  manually,  if  needed,  to  pro¬ 
vide  a  best  fit  to  the  normalized  experimental  penetration  data. 
Results  of  the  parameter  fitting  are  listed  in  the  rightmost  two 
columns  of  Table  3;  here  S  is  multiplied  by  mass  density  to  provide 
a  result  with  the  same  dimensions  as  stress.  Note  that  velocity 
dependent  parameters  are  not  available  for  the  data  provided  in 
Ref.  [5]  because  the  penetration  data  for  the  bare  metallic  backing 
was  not  reported  in  that  study. 

Accuracy,  or  possible  lack  thereof,  of  the  parameter  fits  is 
demonstrated  by  representative  results  in  Fig.  3.  Specifically  shown 
in  each  part  of  Fig.  3  are  a  result  for  each  of  the  six  ceramic  mate¬ 
rials,  encompassing  data  from  all  five  references  in  Table  1  from 
which  both  ^  and  S  could  be  determined.  Velocities  Vo  are  rela¬ 
tively  constant  for  each  distinct  segmented  curve  shown  in 
Fig.  3(a)-(d)  wherein  corresponding  experiments  Vq  was  varied 
over  a  substantial  range,  allowing  for  multiple  such  curves.  In 
Fig.  3(e)  and  (f),  the  range  of  impact  velocities  considered  was  too 
small  to  warrant  construction  of  separate  curves;  instead,  all  data 
are  shown  along  with  linear  fits  that  compare  favorably  despite 
substantial  scatter  among  individual  experimental  data  points. 
Agreement  between  model  and  experiment  is  considered  close 
except  for  the  case  of  titanium  diboride  in  Fig.  3(d)  wherein  the 
model  tends  to  overpredict  penetration  depth  at  the  lowest  impact 
velocity  and  underpredict  depth  at  the  highest  velocity. 

4.  Material  behavior 

As  a  corollary  to  Buckingham's  pi  theorem,  if  a  dimensional 
equation  such  as  (2.1)  relates  a  single  independent  and  dependent 
variable  (n  =  2,  /<  =  1),  then  both  variables  must  have  the  same 
physical  dimensions,  and  their  dimensionless  ratio  must  equal  a 
constant,  i.e., 

y  =/te)=*-^i  =y/X2  =  constant.  (4.1) 

The  present  objective  is  exploration  of  physical  meanings  of 
parameters  ^  and  S  introduced  in  §2  and  fit  to  ballistic  test  data  in 
§3.  Specifically  sought  are  correlations  among  these  parameters 
(which  vary  among  ceramic  materials  and  experimental  data  sets) 
and  measurable  physical  properties  such  as  those  listed  in  Table  2. 
Because  3  and  vary  among  experimental  test  configurations. 


Table  1 

Ballistic  penetration  experiments. 


Experiment 

Projectile 

Lo  [mm] 

Vo  [km/s] 

Backing 

h  [mm] 

Tile  materials 

Rosenberg  &  Yeshurun  [5] 

AP,  steel  core 

12.7,  14.5 

0.92-0.98 

Al  2024-T351 

6-10 

AbOj,  AbO^,  B4C,  SiC,  TiB2 

Hohler  et  al.  [28] 

Rod,  sintered  W 

50,  72.5 

1.2-3.0 

HH  steel 

10-100 

AbO^ 

Sent  et  al.  [29] 

Ogive,  W-C  core 

30 

0.6-1.0 

RHA  steel 

6.7-30 

AbOs ,  AbOs 

Reaugh  et  al.  [30] 

Cylinder,  sintered  W 

25.4 

1.35-2.65 

4340  steel 

6.2-60 

AbOT,  AbO^,  AIN,  B4C,  SiC,  TiB2 

Moynihan  et  al.  [31] 

AP,  steel  core 

35.3 

0.82-0.86 

Al  5083-H131 

1. 3-6.4 

AbO^,  B4C,  SiC 

Savio  et  al.  [23] 

AP,  steel  core 

28.4 

0.60-0.83 

Al  7017 

5.1-9.6 

B4C 
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Table  2 

Representative  ceramic  material  properties. 


Material 

Density 

P  [g/cm^] 

Modulus 

E  [GPa] 

Poisson 
ratio  V 

Tough  Kc 
[MPayfh] 

Compress 
ac  [GPa] 

Bend 
(Tb  [GPa] 

Hard 

Hv  [GPa] 

HEL 

an  [GPa] 

Dyn.  str. 
2r  [GPa] 

Structure  c,a  [nm] 

Inelastic 

deformation 

References 

AI2O3 

3.4 

221 

0.20 

3.0 

2.2 

0.32 

9.1 

5.5 

4.3 

Trigonal  1.3,  0.5 

Slip,  twinning 

[30,36,37] 

AI2O3 

3.8 

373 

0.23 

4.5 

2.6 

0.38 

14.1 

7.6 

5.3 

(same) 

(same) 

[4,10,29] 

AIN 

3.2 

315 

0.24 

2.7 

2.1 

0.35 

11.5 

9.4 

6.0 

Hexagonal  0.5,  0.3 

Slip,  phase  change 

[9,38,39] 

B4C 

2.5 

461 

0.17 

3.1 

2.8 

0.40 

31.4 

16.0 

7.1 

Trigonal  1.2,  0.6 

Amorphization 

[10,36,40] 

SiC 

3.2 

453 

0.16 

5.1 

3.4 

0.40 

27.4 

15.7 

11.4 

Hexagonal  1.5,  0.3 

Slip 

[10,36,41] 

TiB2 

4.5 

521 

0.10 

5.5 

3.0 

0.38 

25.0 

15.0 

13.0 

Hexagonal  0.3,  0.3 

Slip 

[10,36,42] 

Table  3 

Parameter  fits  from  dimensional  analysis. 

Experiment 

Co 

Vr  [km/s] 

r2 

Ceramic 

^  [mm] 

S’  •  p  [GPa] 

Hohler  et  al.  [28]  (Lq  =  72.5  mm)  (Lq  =  50.0  mm) 

-0.67 

1.11 

1  (two  points) 

AbO^ 

120.8 

0.38 

0.57 

3.13 

1  (two  points) 

AbO^ 

62.5 

2.32 

Senf  et  al.  [29] 

-0.44 

0.43 

0.9983 

AbOs 

50.0 

0.019 

AbO^ 

21.4 

0.031 

Reaugh  et  al.  [30] 

0.31 

1.72 

0.9637 

AbOT 

36.3 

-0.81 

AbO^ 

31.8 

-0.15 

AIN 

23.1 

-0.13 

B4C 

25.2 

-0.10 

SiC 

33.0 

1.13 

TiB2 

25.4 

0.72 

Moynihan  et  al.  [31] 

-0.58 

0.45 

0.9882 

AbO^ 

3.78 

-0.12 

B4C 

3.09 

-0.09 

SiC 

3.25 

-0.06 

Savio  et  al.  [23] 

-1.03 

0.28 

0.9997 

B4C 

7.2 

33.0 

Rosenberg  &  Yeshurun  [5] 

AbOj 

9.3,  12.2 

(co,  Vr,  S:  not  available) 

AbO^ 

6.9,  8.1 

(^:  12.7  mm,  14.5  mm  projectile) 

B4C 

6.6,  7.0 

SiC 

6.0,  6.6 

TiB2 

5.7,  6.3 

trends  in  their  values  among  different  materials  can  only  be 
deduced  by  considering  values  for  different  materials  for  fixed 
experimental  configurations,  i.e.,  values  fitted  to  data  from  different 
references  cannot,  in  general,  be  quantitatively  compared. 

Let  {^]t  denote  sets  of  physically  measurable  properties  that  are 
constant  for  a  given  ceramic  tile  material.  Examples  are  listed  in 
Table  2  with  supplementary  references,  though  such  lists  may  be 
incomplete,  i.e.,  other  unlisted  properties  may  also  be  important. 
For  a  fixed  experimental  configuration,  sought  are  relations  of  the 
form 

^  — (4.2) 

Without  further  assumptions,  a  dimensional  analysis  of  (4.2)  is 
fruitless.  As  the  simplest  possible  attempt  at  further  analysis,  as¬ 
sume  that  each  function  depends  most  strongly  on  a  single  material 
property,  to  be  determined.  Then  (4.1 )  yields 


^/f2=C2;  (4.3) 

where  are  ceramic  properties  and  q  are  dimensionless  constants, 
with  i  =  1,2. 

First  consider  f  i,  which  must  have  dimensions  of  length.  Noting 
that  ^  enters  only  the  velocity  independent  part  of  (2.9),  a  logical 
choice  would  be  a  ratio  of  static  properties.  The  best  correlation 
among  possibilities  that  can  be  constructed  from  Table  2  is  ob¬ 
tained  from  ?i  =  r /F,  where  E  is  Young's  elastic  modulus  and  E  is 
surface  energy,  related  to  fracture  toughness  Kc  and  Poisson's  ratio  v 
via 

r  =  /c2(l-.2)/(2£).  (4.4) 

Other  combinations  with  dimensions  of  length  -  e.g.,  F jH,  /F, 
/(tb,  etc.  -  fail  to  yield  trends  or  correlations  with  ^  as  significant 


Fig.  2.  Normalized  residual  penetration  versus  normalized  penetrator  velocity  into  bare  backing  metal  with  linear  fits:  (a)  reference  [30]  (b)  reference  [31]. 
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Fig.  3.  Normalized  residual  penetration  versus  normalized  total  ceramic  tile  thickness,  experimental  data  and  model  results:  (a)  AI2O3 ,  data  from  Ref.  [29]  (b)  AI2O3 ,  data  from  Ref. 
[28]  (c)  AIN,  data  from  Ref.  [30]  (d)  TiB2,  data  from  Ref.  [30]  (e)  SiC,  data  from  Ref.  [31]  (f)  B4C,  data  from  Ref.  [23]. 


as  those  reported  for  T /E  in  what  follows.  A  physical  rationale  for 
correlation  between  E jE  and  ^  is  explained  as  follows.  Recall  from 
(2.9)  and  (2.13)  that  ^  scales  the  effect  of  tile  thickness  h  on 
penetration  depth.  Ratio  FjE  can  be  interpreted  as  a  measure  of 
surface  elastic  energy  to  bulk  elastic  energy.  As  the  ceramic  target 
becomes  thinner  (smaller  h),  the  importance  of  surface  energy 
presumably  becomes  greater,  as  tensile  fracture  associated  with 
bending  becomes  more  likely  [19,22].  Volumetric  elastic  energy 
becomes  relatively  more  important  for  thicker  targets,  as  the  ratio 
of  free  surfaces  to  target  volume  decreases. 

Now  consider  ^2*  which  must  have  dimensions  of  energy  per 
unit  mass  or  velocity  squared.  Noting  that  S  only  enters  the  velocity 
dependent  part  of  (2.9),  a  logical  choice  would  include  dynamic 
material  properties.  The  best  correlation  among  possibilities  that 
can  be  constructed  from  Table  2  is  obtained  from  f  2  =  t/p,  where  r  is 
dynamic  shear  strength  as  defined/reported  in  Ref.  [10]  for 
example,  and  p  is  mass  density.  The  combination  anlp  does  not 
provide  as  satisfactory  a  correlation  with  s,  in  agreement  with  the 
assessment  of  [10]  that  reported  no  correlation  between  the  HEL  of 
the  ceramic  tile  and  residual  penetration  depth.  A  physical  rationale 


for  correlation  between  r/p  and  S  is  explained  as  follows.  Recall 
from  (2.9)  and  (2.13)  that  S  scales  the  effect  of  penetrator  velocity 
Vo  on  penetration  depth.  Ratio  r/p  can  be  interpreted  as  a  measure 
of  frictional  sliding  resistance  or  viscous  stress  supported  by  the 
failed  or  comminuted  ceramic  material  [9,43,44].  Presuming  that 
such  resistance/stress  depends  on  penetrator  velocity  would  sug¬ 
gest  correlation  between  dynamic  shear  strength  and  S.  Normali¬ 
zation  by  mass  density  is  necessary  to  enable  ^2  with  the  same 
dimensional  units  as  S. 

Verification,  or  possible  lack  thereof,  of  the  parameter-property 
correlations  is  considered  in  the  context  of  Table  4,  which  only 
includes  ballistic  results  in  which  multiple  materials  have  been 
tested  for  the  same  experimental  configuration.  First  consider  a 
smaller  value  of  this  parameter  denotes  stronger  resistance  to 
penetration  as  tile  thickness  h  is  increased.  For  each  experimental 
data  set,  ceramics  are  listed  from  top  to  bottom  in  order  of 
increasing  Values  of  FjE  are  listed  side-by-side.  Except  for  the 
three  discrepancies  marked  in  bold  font,  ^  and  F jE  increase  in  the 
same  sequence.  Low  purity  alumina  always  has  a  larger  value  of  ^ 
than  high  purity  alumina,  and  both  aluminas  tend  to  have  larger 
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Table  4 

Parameter-property  correlations:  discrepancies  from  general  trends  in  bold  font. 


Experiment 

Ceramic 

^  [mm] 

lO^r/E  [mm] 

Ceramic 

^  [J/g] 

lO-V/p  [J/g| 

Senf  et  al.  [29] 

AhOt 

21.4 

6.90 

AbOs 

7.87 

0.70 

AbOj 

50.0 

8.85 

AbOs 

5.25 

0.63 

Reaugh  et  al.  [30] 

AIN 

23.1 

3.47 

SiC 

358 

1.78 

B4C 

25.2 

2.19 

TiB2 

160 

1.44 

TiB2 

25.4 

5.52 

B4C 

-39 

1.42 

AbO^ 

31.8 

6.90 

AIN 

-39 

0.94 

SiC 

33.0 

6.17 

Ab03 

-39 

0.70 

AI2O3- 

36.3 

8.85 

AbOs 

-238 

0.63 

Moynihan  et  al.  [31] 

B4C 

3.09 

2.19 

SiC 

-20 

1.78 

SiC 

3.25 

6.17 

Ab03 

-33 

0.70 

AbO^ 

3.78 

6.90 

B4C 

-35 

1.42 

Rosenberg  &  Yeshurun  [5] 

TiB2 

6.0 

5.52 

(^:  average  for  two  penetrator  types) 

SiC 

6.3 

6.17 

(rf:  not  available  from  test  data) 

B4C 

6.8 

2.19 

AbOL 

7.5 

6.90 

AbOs 

10.8 

8.85 

values  than  the  other  ceramics.  Aluminum  nitride  data  is  available 
from  only  one  data  set  [30]  but  corresponds  to  the  lowest  value  of 
for  that  set.  Titanium  diboride  always  has  a  larger  value  than  silicon 
carbide.  Data  for  boron  carbide  are  inconsistent  and  represent  two 
of  the  three  anomalies  from  the  proposed  correlation. 

Now  consider  s,  where  larger  values  denote  improved  pene¬ 
tration  resistance  (decreasing  P/Lq)  with  increasing  velocity  Vq  in 
(2.9).  In  the  rightmost  three  columns  of  Table  4,  ceramics  are 
organized  within  each  data  set  from  top  to  bottom  in  concert  with 
decreasing  s.  Corresponding  values  of  r/p  are  given  as  well, 
demonstrating  increasing  (or  at  least  stationary  values  of)  S  with 
increasing  r/p  except  for  the  lone  data  point  in  bold.  In  two  of  three 
data  sets,  alumina  has  the  lowest  value  (i.e.,  worst  performance), 
and  low  purity  alumina  always  has  a  lower  value  than  high  purity 
alumina.  Silicon  carbide  demonstrates  superior  performance  to  the 
other  ceramics  in  both  data  sets  in  which  it  appears.  The  lone 
anomaly  from  the  correlation  again  pertains  to  boron  carbide, 
which  has  a  lower  value  of  S  but  higher  value  of  r/p  than  alumina 
considered  in  Ref.  [31]. 

Perhaps  unsurprisingly,  three  of  four  deviations  from  proposed 
correlations  arise  for  boron  carbide,  whose  dynamic  behavior  is 
notably  difficult  to  characterize  [45,46].  Some  experimental  con¬ 
figurations  may  lead  to  stress-induced  amorphization  in  the  tile 
while  others  may  not;  [5,30,31]  do  not  provide  microstructural 
characterization  of  the  perforated  tiles  sufficient  to  support  or 
reject  this  conjecture. 

Experimental  data  exist  for  depth  of  penetration  versus  time  for 
confined  semi-infinite  ceramic  targets  impacted  by  high  velocity 
tungsten  long  rods,  specifically  confined  targets  of  alumina  [47], 
aluminum  nitride  [48],  boron  carbide  [49],  and  silicon  carbide  [50]. 
Impact  velocities  were  systematically  varied  from  approximately 
1.5  to  5.0  km/s  in  these  experiments,  with  all  other  factors  apart 
from  the  ceramic  material  essentially  held  fixed.  Though  particular 
Equations  (2.3)-(2.13)  would  not  apply  towards  this  data  (since 
there  is  no  backing  metal),  a  dimensional  analysis  similar  to  that 
invoked  in  the  present  paper  could  be  applied  to  analyze  such  data, 
perhaps  also  incorporating  parameters  similar  to  ^  and  S.  Such  an 
exercise  might  give  further  insight  into  validity  of  the  parameter 
property  correlations  proposed  here  in  §4  of  the  present  paper. 

The  dimensional  analysis  and  resulting  dimensionless  equations 
for  penetration  depth  developed  in  this  paper  provide  a  framework 
for  systematic  future  research.  A  suite  of  additional  ballistic  ex¬ 
periments  is  recommended  that  would  supplement  those  consid¬ 
ered  already,  consisting  of  the  same  penetrator  and  backing 
material  in  each  experiment,  and  varying  the  ceramic  material, 
ceramic  tile  thickness,  and  penetrator  velocity.  Material  properties 


should  be  measured  as  necessary  such  that  uncertainties  in  these 
properties  are  eliminated,  and  then  property-performance  corre¬ 
lations  posited  originally  in  this  work  can  be  more  strongly  vali¬ 
dated,  or  possibly  refined/adjusted  or  refuted.  Numerical 
simulations  as  in  Refs.  [2,9,30]  are  suggested  as  a  means  to  provide 
further  insight  into  possible  correlations  since  physical  properties 
can  be  varied  independently,  and  at  low  cost,  among  simulations  on 
(hypothetical)  ceramic  materials  that  may  not  be  readily  available 
for  ballistic  testing.  Mesoscale  (i.e.,  microstructure  resolution) 
simulations  relating  structure  and  properties  or  property  ratios,  for 
example  as  described  in  Refs.  [51-60],  can  then  be  used  to  suggest 
links  between  material  property  ratios  entering  the  dimensional 
analysis  and  features  of  microstructure  such  as  grain  size,  grain 
boundary  strength,  crystallographic  texture,  etc.  Eor  example,  of 
particular  relevance  to  dimensional  parameter  s  is  previous  work 
in  which  mesoscale  simulations  of  polycrystals  were  used  to 
compute  dynamic  shear  strength  (e.g.,  r)  [55,58].  An  overall  strat¬ 
egy  can  be  written  as 

ballistic  ^  dimensional  ^  material  ^  microstructure 
performance  parameters  properties  (grain  size, 
(P/Lq)  (^,  <^)  (riE.rlp)  orientation...) 

where  arrows  -^denote  links  or  correlations.  Links  labeled  1  and  2 
have  been  formulated  in  the  present  paper;  link  3,  requiring 
consideration  of  microstructure  characterization  data  (experi¬ 
mental  and/or  results  from  aforementioned  multiscale  simulations) 
remains  to  be  addressed  in  more  detail  in  future  work. 

5.  Conclusions 

Principles  of  dimensional  analysis  have  been  applied  towards  a 
study  of  ballistic  penetration  resistance  of  ceramic  materials.  In 
particular,  data  of  study  involve  residual  penetration  depths  into 
thick  metallic  backing  plates  fronted  by  ceramic  tiles  of  alumina 
(low  and  high  purity),  aluminum  nitride,  boron  carbide,  silicon 
carbide,  or  titanium  diboride.  Data  sets  from  six  independent 
experimental  investigations  reported  in  the  literature  have  been 
analyzed. 

Application  of  Buckingham's  pi  theorem  along  with  several 
physical  assumptions  has  led  to  a  dimensionless  penetration  depth 
relation  depending  on  penetrator  velocity,  ceramic  target  thickness, 
and  two  material  parameters  unique  to  the  ceramic  material  when 
other  aspects  of  the  test  configuration  (e.g.,  penetrator  geometry 
and  backing  material)  are  held  fixed  in  an  experimental  series. 
These  parameters  are  a  material  length  scale,  which  is  related  to  the 
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thickness  dependence  of  penetration  resistance  of  the  tiles,  and  an 
energy  per  unit  mass,  related  to  the  velocity  dependence  of  pene¬ 
tration  resistance.  Values  of  the  two  parameters  have  been  deter¬ 
mined  via  fitting  the  dimensionless  relation  to  data  sets  for  the  six 
materials  from  the  experimental  investigations.  Comparison  of 
trends  in  values  among  different  materials  in  a  given  investigation 
with  trends  in  conventionally  measured  ceramic  material  proper¬ 
ties  has  suggested  correlation  between  the  length  scale  and  surface 
(fracture)  energy  divided  by  elastic  modulus  and  correlation  be¬ 
tween  the  energy  per  unit  mass  and  dynamic  shear  strength 
divided  by  initial  mass  density.  Suggestions  for  future  in¬ 
vestigations  involving  numerical  simulations  at  multiple  length 
scales  have  been  outlined  that  would  provide  further  insight  and 
more  definitive  relationships  among  structure,  properties,  and 
performance  in  the  context  of  the  dimensionless  framework. 
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